AD-A078  058  MASSACHUSETTS  INST  OF  TECH  CAMBRIDGE  RALPH  M  PARSONS— ETC  F/6  10/2 
HAVE  POWER  EXTRACTION  BY  FLOATING  BODIES. (U) 

NOV  79  C  C  MEI  r  J  N  NEWMAN 


UNCLASSIFIED 


N0001A-76-C-021* 

NL 


BflC  FILE  con  AD  A 07  8 0  58  Q 


Unclassified _  _ 

SECURITY  CLASSIFICATION  Oc  This  PAGE  M  B  £  1^ 

REPORT  DQCUMENTATio^M  |g  ft  Tp7 

1.  REPORT  NUMBER  1  iRL  GoB?C  AStEsBiTNO.  3.  R£ClP|jfffflkCUTALOG  NUMBER 


I.  REPORT  NUMBER 


14.  TITLE  ( and  Subtitle) 


I  Wave^power  ^extraction  by  floating  bodies 


REPORT  ft  PERIOD  COVERED 

Q  Technical  Impart  , 

— 4v - - - L - — i L - 

PERFORMING  ORG.  REPORT  NUMBER 


[7.  AUTHOR^ 


K.  C.G.^Mei-***  J.N  J  Newman 

3.  PERFORMING  ORGANIZATION  NAME  ANO  AQORESS 

+ 

^  R.M.  Parsons  Laboratory  ,*  *Dept .  of  Civil  Engrg. 
'  Massachusetts  Institute  of  Technology 

Cambridge,  MA  02139  _ _ 


3.  CONTRACT  on  GRANT  NUMBER  (a) 


isJw1, 


014-76-C-0214 

/ 


13.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  »  WORK  UNIT  NUMBERS 


NR  062-228 


II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Office  of  Naval  Research 


- - -  '13.  NUMBER  OF  PAGES 


1 4.  MONITORING  AGENCY  NAME  a  AD ORESSfit  differed Mrom  Controlling  Office)  IS.  SECURITY  CLASS,  (of  thl a  report) 


IS*.  OECLASSI  F|  CATION/ DOWN  GRADING 
SCHEDULE 


|IS.  DISTRIBUTION  STATEMENT  (of  thl  a  Report) 


Material  unclassified;  distribution  unlimited. 


I  17.  DISTRIBUTION  STATEMENT  (of  the  abstract  entered  In  Block  20,  It  different  from  Report) 


DEC  12  1979 


18.  SUPPLEMENTARY  NOTES 


[  19  KEY  WORDS  ( Continue  on  reverse  elde  It  neceeeary  and  Identify  bj  block  number) 


20.  ABSTRACT  (Continue  on  reverse  elde  It  neceeeary  and  Identify  by  block  number) 

'The  linearized  theory  of  water  waves  is  reviewed  in  the  context  of  wave-energy 
absorbers.  A  global  analysis  is  outlined, -leading  to  relative'&imple  expressions 
for  the  maximum  energy  absorption.  Calculations  are  presented  for  the  perform¬ 
ance  of  a  Salter  cam  in  two  dimensions,  and  of  a  Cockerell  raft  in  both  two  and 
three  dimensions.  Similar  rates  of  optimum  energy  absorption  are  anticipated 
for  other  devices  including  the  Kaimei  ship  and  French  air  bag. 


DD  t'™n  1473  EOITION  OF  1  NOV  S3  IS  OBSOLETE 
S/N  0102-LF-014-6601 


" _ Unclassified  y  .  7 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  (HTimi  Data  Bntalad) 


Oii'ut'uaxuM  UN  WAVE  ENERGY  UTILIZATION 


I] 


1. 


INTRODUCTION 


The  conversion  of  wave  energy  from  the  ocean  is  an  important  and 
contemporary  field  of  application  for  the  theoretical  methods  which  have 
been  utilized  in  ship  hydrodynamics  and  other  related  fields.  Since 
there  is  no  extensive  base  of  empirical  design  information  for  wave-power 
devices,  theoretical  methods  are  especially  useful  in  this  field. 

The  basic  theory  for  analyzing  wave-body  interactions  is  well 
developed,  under  the  USual  assumptions  of  small  oscillatory  motions  and 
neglect  of  viscous  effects.  This  approach  has  been  pursued  extensively 
in  the  field  of  ship  hydrodynamics  (Wehausen,  1971;  Newman,  1978)  .  The 
necessary  generalizations  to  analyze  various  wave-energy  devices  require 
trivial  extensions  from  six  degrees  of  rigid-body  motion  to  account  for 
the  pertinent  "higher-Otder  modes"  such  as  the  relative  angular  displace¬ 
ment  at  the  hinge  of  a  Cockerell  raft,  or  the  oscillatory  displacement  of 
the  internal  free  surface  in  a  National  Engineering  Laboratory  device. 

The  classic  theoty  for  wave-body  interactions  is  based  upon  the  use 
of  a  linearized  velocity  potential.  The  overall  problem  can  be  decomposed 
into  linear  components  including  the  incident  wave  system,  the  scattering 
disturbance  of  the  incident  wave  by  the  stationary  body,  and  the  radia¬ 
tion  disturbances  due  to  each  mode  of  body  motion.  For  a  given  body 
geometry  the  scattered  And  radiation  potentials  may  be  derived  by  a 
variety  of  techniques}  generally  involving  numerical  analysis  (Mei, 

1978) . 

In  the  direct  ap^foach  each  of  these  potentials  is  found  and  the 
corresponding  components  of  the  pressure  follow  from  Bernoulli's  equa¬ 
tion.  The  hydrodynamic  forces  are  determined  directly  by  integration 
over  the  wetted  surface  of  the  vessel.  (The  term  "force"  is  used  here  in 
a  generalized  sense,  fcd  include  the  moment  and  any  pertinent  higher-order 
loads.)  Equations  of  HiOtion  follow  for  the  various  degrees  of  freedom 
of  the  vessel,  and  it  is  straightforward  to  solve  for  the  unsteady  mo-  — 
tions  and  for  the  rate  of  energy  absorption. 

An  alternative  gibbai  analysis  can  be  based  in  part  on  the  principle 
of  energy  conservation  and  more  generally  on  the  application  of  Green's 
theorem  to  the  varioud  Separate  components  of  the  velocity  potential. 

This  approach  is  useftlt  particularly  to  determine  the  optimum  power  out¬ 
put  from  a  given  configuration.  In  this  way  it  is  possible  to  determine 
the  maximum  output,  a$t  A  rational  design  goal,  and  to  avoid  the  limita¬ 
tions  of  performance  Mtid  lengthy  experiments  which  result  from  a  purely 
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The  global  analysis  provides  a  relatively  simple  prediction  of  the 
power  absprption,  and  other  characteristics  of  the  device,  including  the 
wave  damping  and  exciting  forces,  in  terms  of  the  far-field  asymptotic 
form  of  the  radiation  potentials.  However,  the  added-mass  force  in  quad¬ 
rature  with  the  body  velocity  must  be  determined  from  a  direct  local  sol¬ 
ution  of  the  radiation  potential.  Thus  a  direct  solution  of  the  radia¬ 
tion  problem  is  required,  for  each  pertinent  degree  of  freedom,  to  pro¬ 
vide  a  complete  description  of  the  relevant  hydrodynamic  characteristics 
of  the  device. 

The  theory  is  simplified  if  the  fluid  motion  is  two-dimensional,' 
implying  a  long  cylindrical  device  with  its  axis  parallel  to  the  wave 
crests.  The  resulting  analysis  has  given  much  inspiration  to  this  field, 
but  practical  considerations  suggest  that  alternative  configurations  with 
smaller  transverse  dimensions  may  be  more  economical,  especially  in 
respect  to  structural  costs  and  mooring  forces. 

This  paper  is  intended  to  review  the  basic  theoretical  approach  to 
wave-power  devices,  and  to  present  a  variety  of  results  for  specific  con¬ 
figurations  which  have  attracted  the  authors’  attention.  The  essential 
elements  of  the  theory  are  outlined  in  Section  2,  with  principal  emphasis 
on  the  rate  of  energy  absorption.  Numerical  results  for  specific  two- 
and  three-dimensional  configurations  are  presented  in  Sections  3  to  5. 

In  Section  6  we  discuss  the  pertinent  hydrodynamic  characteristics  of 
different  types  of  wave-power  devices . 

2.  GLOBAL  ANALYSIS 

With  the  fundamental  assumption  of  small  oscillatory  motions,  the 
incident  wave  spectrum  can  be  decomposed  into  separate  sinusoidal  com-  • 
ponents  of  prescribed  frequency  and  direction.  Therefore,  it  is  relevant 
to  consider  the  absorption  of  a  single  regular  plane  progressive  wave 
system,  the  direction  of  which  is  taken  for  convenience  to  coincide  with 
the  horizontal  x-axis.  The  absorbing  vessel  may  be  oriented  in  an 
arbitrary  manner,  but  to  simplify  our  discussion  only  the  case  of  "head- 
seas"  is  considered. 

Deep  water  is  assumed,  and  the  velocity  potential  of  the  Incident 
waves  is  given  by  the  expression 
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product  is  to  be  taketli  The  parameter  A  denotes  the  wave  amplitude,  g 
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the  gravitational  acceleration,  k  =  w  /g  is  the  wavenumber,  and  z  is  the 
vertical  coordinate. 

The  absorbing  device  is  assumed  to  be  a  floating  or  submerged 
vessel,  which  may  be  rigid  or  flexible.  Its  motions  may  be  defined  by 
six  rigid-body  degrees  of  freedom,  with  complex  velocities  v  ^ ,  and  a 
sequence  of  higher-order  "flexible"  motions  which  can  be  defined  by 
extension  of  the  same  dotation  for  j  >  6.  Since  the  solution  is  linear¬ 
ized,  superposition  of  Separate  components  is  valid  and  the  general 
expression  for  the  velocity  potential  takes  the  form 


<J>  =  (Jij.  +  4>s  +  E 


(2.2) 


Here  $  is  the  scattered  potential,  due  to  the  disturbance  of  the  inci- 
j  s 

dent  waves  by  the  stationary  vessel.  The  remaining  components  4)^  are  due 
to  the  various  degrees  of  freedom  of  its  motion.  In  (2.2)  and  hereafter, 
summation  is  over  all  tlOnzero  modes  of  motion. 

Each  component  of  (2.2)  is  governed  by  Laplace's  equation,  subject 
to  linearized  boundary  Conditions  on  the  free  surface  and  on  the  wetted 
surface  S  of  the  vessel •  With  the  exception  of  the  incident  wave  poten¬ 
tial  <J>  ,  each  component  of  (2.2)  also  satisfies  a  radiation  condition 
X 

of  outgoing  waves  at  lfltge  distance  from  the  vessel,  such  that  for 

2  2 

large  values  of  R  =  (x  +  y  ) , 

=  ie-l7T/4(2TTkRr1/2  Hj(0)eikR  ,  (j  =  s,l,2,...)  (2.3) 


H_.(9)  =  “kjj  (~g^  ”  -^j)  {exp  [kz-ik(x  cos  0  +  y  sin  9)]}dS  (2.4) 
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is  the  So-called  Kochil  function,  proportional  to  the  amplitude  of  the 
far-field  radiated  wavF'S  associated  with  the  potential  . 

From  the  linearized  Bernoulli  equation,  the  pressure  field  of  the 
fluid,  and  thus  the  hycihodynamic  forces  acting  upon  the  vessel,  can  be 
expressed  in  a  similar  fbrm  to  (2.2).  The  total  force  in  the  i-th 
direction  due  to  is  the  exciting  force  due  to  the  incident  wave 
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Similarly,  the  force  due  to  <}>£  is  expressed  in  the  form 


<t>  -H-  <?S  =  B.  .  -  iu>  M. . 

j  ij  ij 


(2.6) 


where  the  real  coefficients  B,  .  and  M.  .  are  known  as  the  damping  and 

ij  ij 

added  mass,  respectively.  The  damping  coefficients  are  especially  import¬ 
ant  in  the  context  of  wave  power  •/  as  the  components  of  the  hydrodynamic 
force  in  phase  with  the  velocity  of  the  vessel;  these  are  responsible 
for  the  radiated  wave  energy  due  to  the  vessel's  own  motions  in  calm 
water. 

A  direct  analysis  of  wave  power  absorption  requires  solutions  for 
the  scattered  potential  and  each  of  the  "forced  motion"  potentials  in 
calm  water,  following  which  the  pressure  force  acting  in  each  case  can  be 
found  from  integration.  It  is  clear  however  from  the  "principle  of 
virtual  work"  that  the  power  absorbed  will  depend  only  on  the  hydrodynam¬ 
ic  force  coefficients  (2.5)  and  (2.6),  and  the  mechanical  device  restrain¬ 
ing  the  body  motions. 

Instead  of  solving  all  the  boundary  value  problems,  much  Information 
may  be  obtained  by  using  certain  global  relations  which  follow  mathemati¬ 
cally  from  Green's  theorem.  These  are  examined  in  a  systematic  manner 
by  Newman  (1976) .  The  most  important  results,  in  the  context  of  wave- 
power  absorption,  are  relations  for  the  exciting  force,  damping  coeffi¬ 
cients,  and  absorbed  power,  in  terms  of  the  Kochin  functions,  as 
summarized  below: 


The  Haskind  relations  express  the  exciting  force  in  the  form 

•-Tv,) 


(2.7) 


Here  is  the  complex  amplitude  of  the  exciting  force  in  waves,  of 
unit  amplitude,  with  j  =  1,2,3  for  each  component  of  the  force,  j  «  4,5,6 
for  the  moment.  Higher-order  mode  shapes  and  exciting  forces  can  be 
defined  similarly.  The  principal  advantage  of  (2.7)  is  to  permit  the 
determination  of  the  exciting  force  in  waves  without  solution  of  the 
boundary  value  problem  for  the  scattered  potential  <J>g. 

The  damping  coefficients  are  related  more  directly  to  the  outgoing 
radiated  waves  by  the  principle  of  energy  conservation.  Mathematically 
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this  relationship  tak^s  the  form 
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(2.8) 


where  B_^  denotes  the  damping  coefficient  for  the  force  or  moment  in  the 
direction  i,  due  to  a  body  motion  in  mode  j.  The  matrices  .  and  M^. 
are  symmetric. 

For  any  combination  of  body  modes  such  that  the  phase  of  the  motion 
is  constant  over  the  body  surface,  in  each  mode,  with  complex  amplitude 
Vj ,  the  rate  of  energy  absorbed  by  the  body  can  be  expressed  in  the  form 


n  -  - j*  w  '?>>>  -  n  n*  vh*- 


(2.9) 


The  last  expression  cat!  be  interpreted  physically  by  combining  it  with 
(2.7)  and  (2.8)  to  give 


f  ’  1  A  **tE  V  "  2  vi  XJ  BiJ 


(2.10) 


This  expression  is  equivalent  to  equation  (4.1)  of  Evans  (1979),  where  a 
matrix  notation  is  used  and  the  derivation  is  based  on  an  analysis  of  the 
work  done  in  opposition  to  the  hydrodynamic  forces  on  the  body.  Evans' 
derivation  is  for  a  system  of  N  bodies,  each  having  a  single  degree  of 
freedom,  but  there  is  itt  principle  no  difference  between  that  case  and 
the  present  one,  of  a  single  body  with  several  degrees  of  freedom. 

Indeed,  the  generalization  to,  say,  N  bodies  with  M  degrees  of  freedom 
simply  requires  a  suitable  set  of  additional  subscripts,  with  no  other 
change  to  the  analysis. 

The  above  equatioris  are  valid  in  general  for  a  three-dimensional 
vessel.  In  the  special  case  of  two-dimensional  motions,  resulting  from 
the  interaction  of  a  lh^g  cylindrical  vessel  with  its  axis  parallel  to 
the  wave  crests,  equations  (2.1),  (2.2)  and  (2.8)  remain  valid,  but 
the  intermediate  equations  are  modified  as  noted  by  Newman  (1976)  . 

The  equivalent  expressions  (2.9)  and  (2.10)  can  be  used  to  determine 
the  performance  of  a  wMfe-absorbing  device  simply  in  terms  of  its 
respective  Kochin  functions,  or  in  terms  of  its  damping  coefficients  and 
exciting  forces.  The  feetond  term  in  (2.9)  or  (2.10)  is  a  negative 
quantity,  associated  pH^Hically  with  outgoing  energy  flux  due  to  the 
vessel’s  own  wave  genePSHon.  The  first  term  in  these  expressions  may 
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be  of  either  sign,  depending  on  the  relative  phase  between  the  velocity 
component  and  corresponding  Kochin  function  or  exciting  force. 

Optimum  performance  follows  by  choosing  the  velocity  components  in  a  suit¬ 
able  manner. 

If  there  is  only  one  degree  of  freedom,  the  energy  absorption  is 
maximized  by  choosing  the  phase  of  this  velocity  to  be  the  same  as  tne 
exciting  force,  and  the  resulting  quadratic  function  attains  a  maximum 
value  when 
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(2.11) 


It  can  be  shown  (cf.  Mei,  1976)  that  this  situation  is  realized  if  the 
vessel  is  restrained  by  a  linear  device  tuned  for  resonance  in  such  a 
manner  that  the  inertial  and  restoring  terms  cancel,  and  with  a  mechan¬ 
ical  damping  precisely  equal  to  the  hydrodynamic  damping. 

Corresponding  to  (2.11),  the  maximum  power  absorption  is  given  by 
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(2.12) 


In  this  form  it  is  apparent  that  an  efficient  wave  absorbing  device  is 
one  which  radiates  waves  primarily  in  the  direction  opposite  to  the  wave 
incidence  or,  equivalently,  a  vessel  which  experiences  a  relatively  large 
exciting  force  in  waves  from  this  particular  direction.  In  other  words, 
a  geometrically  efficient  three-dimensional  wave  absorber  should  be  such 
that  it  can  generate  highly  focused  waves  by  its  own  motions. 

The  two-dimensional  analog  of  (2.12)  is  noted  by  Mei  (1976)  and 
Evans  (1976) ,  and  is  essentially  the  mathematical  basis  for  the  superior 
performance  of  unidirectional  wavemakers  as  wave  absorbing  devices.  In’ 
particular  the  optimum  efficiency  which  is  attainable  with  one  degree  of 
freedom  is 


(i  +  |a+/aJ2)-1 


(2.13) 


where  A+  and  A  are  the  radiated  wave  amplitudes  towards  x  ^  -K°  and 

respectively,  the  incident  wave  being  from  x  ^  -°°.  In  this  context  the 

efficiency  is  defined  as  the  ratio  of  the  power  absorbed  by  the  device 

12  2 

and  the  rate  of  energy  flux  ^  pg  A  /u  per  unit  width  of  the  incident 
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wave  system. 

The  results  for  a  6ingle  degree  of  freedom  can  be  extended  to  a  pair 
of  modes  of  motion  which  are  even  and  odd  with  respect  to  the  x-axis,  for 
example  heave  and  surge  in  the  case  of  rigid-bodv  motions.  The  corre¬ 
sponding  Kochin  functions  are  even  and  odd  in  6,  and  these  do  not  inter¬ 
act  in  (2.9).  Equivalently,  in  terms  of  (2.10),  there  is  no  cross¬ 
coupling  between  such  modes  and  B„  =  0  for  i  -  j  .  It  follows  that  the 
total  energy  absorption  is  the  linear  superposition  of  that  in  each 
separate  mode. 


The  three-dimensional  results  may  be  illustrated  by  considering  an 

axisymmetric  vessel.  Heave  motions  are  independent  of  0,  and  from  the 

first  form  of  (2.12)  it  follows  that  the  maximum  energy  absorption  is 
2  2 

pg  A  /4iok.  Dividing  by  the  rate  of  energy  flux  per  unit  width  of  wave 
crest  in  the  incident  waves  gives  a  "capture  width"  W  =  1/k,  or 
approximately  one  sixth  of  the  wavelength.  Alternatively,  surge  or  pitch 
of  the  same  vessel  gives  a  Kochin  function  proportional  to  cos  0,  and 
W  =  2/k.  Superposition  of  these  two  modes  gives  a  maximum  capture  width 
from  rigid-body  motions  of  3/k.  Larger  amounts  of  energy  can  be  captured J 
in  principle,  by  highet-order  modes  proportional  to  cos  n0 .  In  all  cases 
the  optimum  combination  of  an  even  and  odd  mode  shape  is  related  in 
precisely  the  manner  fcd  cancel  the  radiated  waves  downstream,  (0  =  0) 
and  reinforce  the  sepatate  components  upstream  (0  =  tt)  .  . 

Optimum  performance  of  a  wave-absorbing  device  requires  that  each 
mode  is  optimized  in  tHflgnitude  and  phase.  Equivalently,  each  mode  must 
be  restrained  by  a  fotfce  of  optimum  magnitude  and  phase.  Thus  the  re¬ 
straint  which  extracts  the  output  power  must  have  an  optimum  complex 
impedance  including  irieftial  and  dissipative  components.  This  optimum 
impedance  is  a  function  of  the  incident  wave  frequency  or  wavelength.  A 
simple  mechanical  device  with  constant  inertia  cannot  satisfy  this  re¬ 
quirement  for  all  wave ' fcngths .  The  global  analysis  can  be  used  in  non¬ 
optimum  conditions  to  express  the  rate  of  energy  absorption  to  the  body 
motions,  but  the  latter  must  be  computed  by  direct  solution  of  the 


equations  of  motion. 


3.  SALTER'S  CAM 


Salter's  cam  is  tta  first  device  which  aroused  serious  interest  in 
wave  power  by  its  remiffcsbly  high  efficiency  unimagined  and  unattained 
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pressure.  This  implies  a  very  large  force;  e.g.,  for  a  =  10  m,  A  =  3m, 
PgaA  =  30  Tn/m.  It  in  interesting  that  for  extraction  rates  below 
optimum  the  peak  force  and  the  peak  efficiency  do  not  occur  at  precisely 
the  same  frequency.  Mynett  et  al.  (1979)  also  performed  calculations 
for  various  submergence  and  water  depths,  and  found  the  latter  to  have 
relatively  little  effect. 

In  reality  a  rigid  shaft  would  demand  a  massive  support  which  of 
course  implies  a  high  Cost.  A  flexible  support  would  normally  lead  to 
a  loss  in  efficiency.  Calculations  based  on  a  simple  model  in  which  the 
cam  may  slide  on  a  vertical  sleeve  which  may  move  freely  horizontally 
show  that  the  efficiency  is  usually  reduced  drastically,  with  similar 
reduction  in  wave  forces.  This  is  illustrated  in  Figure  3.  In  prin¬ 
ciple,  mechanisms  carl  be  added  in  the  cam  so  that  energy  may  be  extracted 
from  all  degrees  of  freedom  to  increase  the  efficiency,  but  intricate 
designs  would  be  necessary  which  raise  new  technological  and  economical 
considerations.  Once  the  details  of  the  mechanisms  are  specified, 
computations  similar  tt)  those  given  here  are  also  possible. 

If  the  wind-generated  wave  is  assumed  to  be  a  stationary  random 
process  the  averaged  power  flux  per  unit  length  of  the  crest  may  be 
related  to  the  energy  Spectrum  by 


<P  >  =  Pg 
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(3.1) 


while  the  power  extracted  by  the  rolling  cam  is 


,  »  °2 
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where  a^to)  is  the  complex  frequency  response  of  the  roll  amplitude. 

The  mean  squares  of  thfe  roll  amplitude  and  of  the  hydrodynamic  force  on 
the  shaft  can  also  be  Obtained  from  the  incident  wave  spectrum  and  the 
frequency  response: 


'a2>  =  j°°  SjXto)  |u j 1 2dw 

0 
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Serman  and  Mei  (1979)  have  computed  these  mean  squares  based  on  the  ^ 

frequency  responses  o£  Mynett  et  al.  (1979).  „  „  -vrulobla  to  j0®*.  . 
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As  wave  power  devices  will  likely  be  installed  near  the  coastline, 

a  representative  JONSVJAP  spectrum  is  chosen  for  S^..  The  spectrum  is 

characterized  by  the  wind  speed  U  and  the  fetch  length  X  from  which 

various  parameters,  in  particular  the  peak  frequency  and  root  mean 

square  of  the  incident  wave  height,  may  be  determined  by  empirical 

formulas  (see  Hasselmann  et  al.,  1973).  From  JONSWAP  spectrum  the 

increase  of  root-mean- square  wave  height  and  peak  period  with  both  U  and 

X  may  be  calculated.  The  following  numerical  results  for  a  cam  with 

radius  a  =  5  m  are  due  to  Serman  and  are  similar  to  those  given  in 

Serman  and  Mei  (for  a  =  3m).  The  extracted  power  in  kw/m  of  course 

increases  with  both  U  and  X  as  shown  in  Figure  4.  When  U  and  X  increase, 

the  peak  frequency  decreases,  hence  the  efficiency  also  decreases  (Figure 

5); this  is  ct nsistent  with  the  e  curve  in  Figure  1.  Root-mean-square 

opt 

horizontal  and  vertical  forces  are  shown  in  Figures  6  and  7.  The  root- 
mean^square  of  the  angular  displacement  is  shown  in  Figure  8.  For 
example,  let  X  =  200  Km,  U  =  3  ^  14  m/s ec,  <P>  =  2  'v,  46  kw/m,  F^  = 

1  10  Tn/m.  It  should  be  remarked  that  for  these  ranges  of  X  and  U  the 
wave  height  is  2  a,  4  m  so  that  linear  approximation  is  still  useable. 


Serman  and  Mei  (1979)  also  studied  the  performance  in  a  Pierson- 
Moskowitz  spectrum  of  a  cam  whose  shaft  can  move  passively  in  both  the 
vertical  and  horizontal  directions.  Details  are  omitted  here. 


4.  TWO-DIMENSIONAL  FLOATING  RAFTS 

* 

The  elongated  and  shallow  rafts  of  Cockerell  ,  which  extract  energy 
from  the  relative  rotation  about  the  hinges,  have  the  obvious  advantage 
of  smaller  mooring  forces.  It  was  therefore  felt  desirable  to  provide  a 
theoretical  model  for  its  performance.  For  analytical  simplicity  Haren 
and  Mei  (1979)  studied  the  case  of  two-dimensional  rafts  on  a  sea  of 
small  depth.  Because  they  were  also  interested  in  optimum  raft  lengths, 
computations  based  on  the  fully  two-dimensional  theory  for  arbitrary 
wavelength  would  be  quite  tedious.  Therefore  a  linearized  shallow  water 
theory  was  used.  It  is  interesting,  however,  that  for  a  sample  geometry 
the  difference  between  these  two  theories  was  also  found  later  to  be 
not  large  for  a  wide  range  of  frequencies.  Again  according  to  (2.12) 

Essentially  the  same  idea  has  been  proposed  by  Glen  Hagen  of  Louisiana, 
U.S.A.,  in  a  patent  application  (No.  4077213,  Official  Gazette,  U.S. 
Patent  and  Trademark  Office,  March  7,  1978,  filed  Feb.  13,  1976).  The 
elongated  feature  was  stressed  more  by  Hagen  than  by  Cockerell. 


optimum  efficiency  is  attained  if  most  of  the  radiated  waves  are  towards 
x  %  -oo  and  very  little  towards  x  'u  +°°.  This  is  possible  if  the  rafts  / 
on  the  transmission  side  are  much  longer  than  the  rafts  on  the  incidenca\ 
side.  For  the  simplest  case  of  a  short  raft  hinged  onto  a  fixed  semi-  I 
infinite  raft  (which  his  the  same  effect  as  a  fixed  wall),  it  can  be 
easily  shown  that  10dX  efficiency  is  possible  if  the  short  raft  is  res-  U 
onated  and  the  extraction  rate  equals  the  roll  damping  rate  of  the  short 
raft.  With  only  a  buoyancy  restoring  force  the  corresponding  raft  length 
is  40%  of  the  wavelength.  The  efficiency  curve  is  again  quite  broad- 
banded.  For  a  finitd  total  length  without  rigid  external  constraint, 
some  efficiency  is  eApBcted  to  be  lost.  Experimenting  with  three,  four 
and  five  rafts,  HareU  #nd  Mei  sought  the  optimal  design  under  various 
criteria.  The  optimising  parameters  are  raft  lengths  and  extraction 
rates  at  the  hinges.  The  criteria  include  the  best  efficiency  at  a 
fixed  total  length,  while  the  incident  waves  may  be  regular  or  have  a 
smooth  spectrum  with  one  or  two  peaks. 

For  example,  for  a  depth  of  9.15  m,  a  regular  wave  of  period  T  = 

P 

5.52  sec,  wavelength  ■*  52  m,  total  length  =  137  m,  the  best  efficiency 

is  achieved  when  the  three  rafts  are  16.8,  52,  68.3  m  long  and  the 

extraction  rates  at  the  first  two  hinges  are  nearly  equal  to  5380 
2 

kg/m  .  The  total  amount  of  energy  extracted  at  the  first  hinge  exceeds 
95%  of  the  total.  TH*  efficiency  for  other  frequencies  is  given  in 
Figure  9.  It  may  be  Concluded  that  the  relative  notion  about  the  second 
hinge  is  quite  small  90  that  the  system  behaves  like  a  two-raft  system. 
Since  the  entire  syst£ft  is  floating,  the  forces  on  the  hinges  are  not 
large;  a  typical  vertical  force  is: 

F  'u  0.6  pgAH  =  Tn/m  at  resonance,  for  A  =  3  m,  I 

z 

which  is  somewhat  led!  than  that  experienced  by  a  Salter  cam  on  a  fixed 
axis.  In  practice,  ifc  may  be  easier  to  have  a  strong  hinge  between  twp 
moving  rafts  than  to  hftve  a  strong  anchor  to  the  sea  bottom. 

Finally  in  many  hbastal  waters  the  presence  of  swell  introduces  two 
spectral  peaks.  The  *-^sulrs  of  optimization  with  some  consideration  of 
cost  is  rather  interesting.  It  may  be  shown  that  the  piofit  index  is 
<eff>  -  QL^,  where  Q  tS  proportional  to  the  unit  cost  of  the  raft  and 
inversely  proportional  to  the  price  of  electricity  and  the  local  wave 
energy.  Sample  calculations  for  the  Northwest  Pacific  coast  of  the  U.S. 
show  that  the  efficieHhy  curve  follows  the  short  wave  spectral  peak  as 
shown  in  Figure  10,  sWK^esting  that  it  is  cost-effective  to  disregard 
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the  long  wave  peak.  Naturally  a  larger  Q  can  lead  to  negative  profit. 


5.  ELONGATED  SLENDER  VESSELS 

The  final  numerical  example  to  be  presented  is  that  of  a  three-dim¬ 
ensional  Cockerell  raft  which  is  elongated  in  the  direction  perpendicular 
to  the  wave  crests.  More  details  are  given  by  Newman  (1979).  A  similar 
analysis  and  conclusions  should  apply  to  other  devices  including  the 
Kaimei  ship  and  French  air  bag. 

We  consider  a  body  which  is  elongated  in  the  x-direction,  and  sub¬ 
ject  to  a  vertical  displacement  V(x)  along  its  depth.  This  motion  may 
be  expressed  conveniently  in  terms  of  a  sequence  of  modal  functions,  in 
the  form 


V(x)  =  Z  Vj  f  (x) 


C5.1J 


where  v^  is  a  complex  constant 


For  a  slender  body  with  small  beam  and 
draft,  compared  to  the  wavelength,  the  Kochin  function  (2.4)  can  be 


approximated  in  a  consistent  manner  by  neglecting  the  second  term  in 
* 

parenthesis.  For  each  of  the  above  modes,  it  follows  that 


V9) 


-kf 


1/2 


-112 


fj(x)  exp  (ikx  cos  0)  b(x)dx 


(5.2) 


where  Si  is  the  body  length  and  b(x)  is  the  local  waterplane  beam. 


For  each  mode  j,  the  maximum  rate  of  energy  absorption  can  be  deter¬ 
mined  by  substituting  (5.2)  in  (2.12).  In  the  following  examples  the 
local  beam  b  is  assumed  constant,  corresponding  to  a  vessel  with  rec¬ 
tangular  waterplane.  The  rate  of  energy  absorption  is  expressed  as  the 
product  of  the  rate  of  energy  flux,  per  unit  of  incident  wave  front,  and 
the  capture  width  W  of  the  vessel. 

From  the  computational  standpoint,  the  Legendre  polynomial  P^(x)  is. 

a  convenient  choice  of  mode  shape.  The  Kochin  functions  (5.2)  are  given 

in  terms  of  spherical  Bessel  functions.  Moreover,  the  modes  j  =  0,1 

correspond  respectively  to  rigid-body  motions  in  heave  and  pitch,  while 
_ 

This  approximation  is  essentially  identical  to  that  employed  by  Evans 
(1979)  to  study  the  interaction  between  several  adjacent  vessels,  under 
the  assumption  that  each  vessel  operates  independently  without  mutual 
hydrodynamic  interactions.  For  that  reason  Evans’  conclusions  are 
closely  related  to  those  found  here. 


vyn 


higher-order  inodes  correspond  to  flexural  motions  of  polynomial  form. 

The  capture  width  W  corresponding  to  each  of  the  first  six  Legendre 
polynomials  is  shown  by  the  corresponding  solid  curve  in  Figure  11.  A 
common  feature  is  that  for  long  wavelengths,  W  is  proportional  to  the 
wavelength.  This  follows  from  the  results  noted  at  the  end  of  Section  2 
and  from  the  fact  that,  if  X/8,  »  1,  the  vessel  geometry  is  not  signifi¬ 
cant.  However,  the  magnitude  and  phase  of  each  node  is  optimized  in 
these  results,  on  5  basis  of  (2.11),  and  for  long  wavelengths  the 

corresponding  values  of  (2.11)  are  unrealistically  large.  Indeed,  it 

i+2 

is  readily  shown  in  this  case  that  v^  =  0(AJ  )  for  large  values  of  A. 

From  the  practical  standpoint  it  is  appropriate  to  limit  the  magni¬ 
tude  of  the  body  motions,  to  be  of  the  same  order  as  the  incident  wave 
amplitude.  Two  pragmatic  limits  are  shown  by  dashed  lines  in  Figure  11. 
The  relevant  parameter  here  is  the  product  of  v^  and  the  beam-length 
ratio  of  the  vessel.  Representative  values  for  this  product  have  been 
chosen  to  correspond  with  typical  values  of  the  beam-length  ratio  for 
a  slender  wave-absorbet  of  ship-like  proportions. 

The  extraction  of  wave  power  from  the  rigid-body  motions  in  pitch 
and  heave  requires  an  Bkternal  force  or  moment  acting  on  the  vessel,  and 
therefore  a  rigid  foundation  or  mooring  system.  The  higher-order  modes 
of  motion  are  attractive  in  so  far  as  power  can  be  extracted  from  the 
relative  motion  between  different  elements  of  the  device’,  which  is 
moored  only  against  the  second-order  mean  drift  force  of  the  waves.  For 
this  reason  it  is  appropriate  to  impose  the  condition  that  the  vessel  is 
free  to  respond  to  the!  Waves  in  heave  and  pitch,  without  an  external 
first-order  restraint. 

As  a  consequence  df  the  slender-body  approximation,  the  latter 
condition  is  satisfied  by  the  higher-order  modes  n  =  2,3,...,  and  more 
generally  by  any  mode  shape  which  is  orthogonal  to  the  Legendre  poly¬ 
nomials  Pq  =  1  and  Pj  -  x. 

Trigonometric  modfe  shapes  tuned  to  the  incident  wavenumber  can  be 
expected  to  give  betthi:  performance  than  a  simple  polynomial,  expecially 
for  shorter  wavelengths.  A  suitable  pair  of  even  and  odd  modes  ortho¬ 
gonal  to  heave  and  pi^th  are  given  by  the  expressions 


f  =  cos(kx)  -  j^fk.1/2), 

fQ  =  sin(kx)  -  (d#/Jl)  j1(kJL/2)  , 
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Here  jg  ^  are  spherical  Bessel  functions  of  the  first  kind.  The  corre¬ 
sponding  capture  widths  are  shown  in  Figure  12;. .  These  results  are 
similar  to  the  Legendre  function  modes  n  -  2,3  for  long  wavelengths,  but 
for  short  wavelengths  the  improved  performance  is  apparent. 

As  an  alternative  to  these  continuous  mode  shapes  we  consider  an 
articulated  raft  with  rigid  elements,  joined  by  two  symmetrically  placed 
hinges.  TWo  modes  of  motion  are  defined,  orthogonal  to  heave  and  pitch, 
as  depicted  in  Figure  13.  The  corresponding  values  of  the  capture  width 
are  shown  in  the  same  figure.  Similar  calculations  indicate  that  the 
hinge  position  is  not  an  important  parameter,  provided  only  that  the 
hinges  are  not  close  to  the  ends  of  the  raft.  A  raft  with  three  elements 
of  equal  length  appears  to  be  most  practical.  It  follows  also  that  a 
raft  with  only  one  hinge  at  the  center  will  perform  with  approximately 
the  same  performance  as' for  the  even  mode  shown  in  Figure  12,  but  without 
the  added  degree  of  freedom  permitted  by  the  corresponding  odd  mode. 

(When  the  two  hinges  approach  the  center,  the  odd  mode  shown  in  this 
figure  tends  to  a  discontinuous  shearing  motion  at  the  center.) 

Since  the  even  and  odd  modes  can  be  superposed,  the  total  capture 
width  in  the  above  examples  may  approach  a  peak  value  close  to  the  total 
length  i  of  the  vessel.  Since  the  continuous  mode  shapes  shown  in 
Figures  11  and  12  attain  similar  capture  widths,  it  is  likely  that  the 
articulated  raft  with  two  hinges  is  nearly  optimum  and  no  significant 
benefit  will  result  from  the  number  of  hinges  beyond  this  nunfcer. 

With  optimum  tuning  the  even  and  odd  modes  of  body  motion  are  in 
quadrature.  These  combine  to  give  a  wavelike  motion  along  the  vessel 
with  the  phase  velocity  of  the  incident  wave.  The  amplitude  of  this 
motion  is  symmetrical  fore-and-aft,  in  agreement  with  the  analogous 
finding  of  Evans  (1979).  This  result  differs  from  the  conclusion  in 
Section  4  based  on  a  two-dimensional  mode,  where  the  optimum  motion  is  . 
concentrated  primarily  at  the  end  facing  the  incident  waves .  This  dis¬ 
tinction  can  be  attributed  to  three-dimensional  diffraction,  with  the 
Incident  wave  energy  focused  upon  the  after  end  of  the  body  from  the 
sides. 

The  performance *ln  a  spectrum  can  be  computed  as  in  (3.1)  -  (3.2). 
As  an  example  we  show  in  Figure  13  the  total  power  extracted  in  a  Pier- 
son-Moscowitz  spectrum  by  a  raft  with  two  hinges,  as  a  function  of  the 
overall  raft  length. 
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The  calculations  df  this  section  are  based  on  the  "ordinary"  slender 
body  approximation  which  assumes  the  beam  and  draft  small  compared  not 
only  to  the  length  of  the  vessel,  but  also  to  the  wavelength.  A  "uni¬ 
fied"  approximation  ha9  been  developed  to  remove  the  latter  restriction, 
as  described  by  Nays  (1978)  and  Newman  (1978) .  Recent  computations  with 
this  more  refined  approach  indicate  capture  widths  which  are  slightly 
less  than  those  shown  in  Figure  12,  but  this  reduction  is  at  most  10%. 
Further  work  to  improve*  upon  the  slender-body  approximation  will  be  re¬ 
ported  in  future  publications . 


6.  DISCUSSION  AND  CONCLUSIONS 


Devices  for  absorbing  wave  energy  can  be  categorized  from  the  hydro- 
dynamic  viewpoint  accdfding  to  the  following  list: 

Cl)  Floating  or  submerged 

(2)  Rigid  or  flexible  elements 

(3)  Stationary  Ot  freely  floating 

(4)  Two-  or  thred^dimensional  configuration 

(5)  Large  devicd  Or  point  absorber 

The  first  two  characteristics  in  this  list  are  of  relatively  little  im¬ 
portance  from  the  standpoint  of  wave-power  absorption.  This  is  suggested 
by  the  global  analysis  in  Section  2,  from  which  it  follows  that  the  op¬ 
timum  power  absorption  depends  essentially  on  the  ability  of  the  vessel 
to  radiate  waves  in  a  Preferred  direction.  As  in  the  analogous  case  of 
laboratory  wavemakers*  this  property  can  be  achieved  with  more-or-less 
equal  success  by  a  variety  of  devices  which  are  floating  or  submerged, 
consisting  of  either  figid  elements,  pneumatic  chambers,  or  other  means 
for  forcing  oscillatory  motion  of  the  adjacent  fluid. 

A  comparison  of  figures  1  and  10  illustrates  the  similarity  in 
optimum  performance  fH{-  two  different  types  of  device,  the  Salter  cam  { 
and  the  two-dimensional  Cockerell  raft.  Similar  results  may  be  antici¬ 
pated  for  all  two-dim»Hsional  devices.  Indeed,  any  two-dimensional  ves¬ 
sel  with  two  or  more  *>h: table  degrees  of  freedom  and  an  optimal  control 
system  can  operate  witii  100%  efficiency! 


More  important  ill  fche  question  of  whether  the  device  is  stationary 
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or  freely  floating.  In  this  context  "stationary"  implies  a  fixed  refer¬ 
ence  from  which  to  exert  an  optimum  restraining  force  on  the  device. 

Most  wave  absorbing  devices  have  been  designed  on  this  basis.  Generally, 
if  they  are  freely  floating,  or  fixed  in  mean  position  by  a  "slack"  moor¬ 
ing  which  exerts  no  first-order  oscillatory  force,  the  resulting  perfor¬ 
mance  is  degraded.  Figure  3  illustrates  this  effect  for  the  case  of  a 
Salter  cam.  Large  structures  which  are  mounted  on  the  ocean  bottom  in  a 
fixed  manner  should  be  included  in  the  same  category  as  stationary  ves¬ 
sels.  By  comparison,  devices  which  extract  power  through  the  relative 
motion  between  adjacent  elements  can  be  designed  for  efficient  perform¬ 
ance  with  a  slack  mooring  system. 

The  fourth  characteristic  listed  above  is  intended  to  place  in  a 
separate  category  those  devices  which  are  based  on  a  two-dimensional 
analysis,  and  which  therefore  must'be  elongated  in  the  direction  parallel 
to  the  wave  crests.  This  two-dimensional  configuration  offers  the  obvi¬ 
ous  advantage  of  absorbing  power  uniformly  along  a  broad  expanse  of  the 
incident  wave  front.  In  principle,  100%  efficiency  can  be  attained  in 
this  manner,  and  the  theoretical  analysis  associated  with  these  devices 
is  relatively  simple.  On  the  other  hand,  two-dimensional  configurations 
must  expose  a  large  projected  area  to  the  incident  waves,  and  therefore 
experience  proportionally  large  mooring  forces. 

Three-dimensional  devices  interact  with  the  surrounding  wave  field 
in  a  more  complicated  manner.  Suitably  designed,  such  devices  have  the 
potential  ability  to  absorb  the  wave  energy  over  a  capture  width  substan¬ 
tially  larger  than  the  transverse  projection  of  the  device  itself.  This 
desirable  attribute  follows  from  the  nature  of  three-dimensional  wave 
diffraction,  which  is  enhanced  by  the  optimum  interactions  of  different 
elements  of  the  vessel.  Analogous  effects  result  from  the  mutual  inter¬ 
action  of  several  small  vessels,  in  the  manner  described  by  Evans  (1979) . 
Similar  considerations  are  well  known  in  the  design  of  radio  antennas 
and  phased  sonar  arrays. 

The  last  of  the  hydrodynamic  characteristics  in  the  list  above  is 
the  size  of  the  vessel  in  relation  to  the  wavelength.  Vessels  which  are 
long  can  diffract  the.  incoming  waves  in  the  manner  described  above, 
while  those  which  are  wide  can  operate  in  a  quasi-two-dimensional 
manner.  Intermediate  between  these  are  axisymmetric  vessels,  of  inter¬ 
est  primarily  in  the  context  of  small  "point-absorbers”  with  horizontal 
dimensions  much  less  than  the  wavelength.  The  economy  and  operational 


simplicity  of  such  devices  has  been  stressed  especially  in  the  work  of 
Budal  and  Falnes.  Point  absorbers  take  advantage  of  three-dimensional 
diffraction  to  absorb  the  incident  energy  over  a  maximum  capture  width  of 
half  a  wavelength,  as  noted  in  Section  2.  (For  heaving  motion  only,  the 
maximum  capture  width  is  about  one  sixth  of  a  wavelength.)  Large  ampli¬ 
tudes  of  motion  are  necessary  to  approach  these  maximum  capture  widths, 
and  a  stationary  reference  frame  is  required  for  the  restraining  force. 

The  most  important  hydrodynamic  characteristics  of  a  wave  absorbing 
device  are  its  size  and  configuration  in  relation  to  the  incident  wave 
field,  and  its  requirement  for  a  stationary  reference  position  or  alter¬ 
natively  a  slack  mooring.  There  are  obvious  practical  and  theoretical 
advantages  associated  with  an  elongated  device,  such  as  a  slender  Cocker¬ 
ell  raft.  The  same  considerations  apply  to  other  devices  with  a  similar 
configuration,  notably  the  Kaimei  ship  and  the  French  air  bag.  With 
suitable  optimization  of  their  design,  these  elongated  vessels  are  capabl* 
of  performing  at  similar  levels,  with  capture  widths  comparable  to  their 
lengths. 

The  computations  in  Section  5  are  intended  to  illustrate  the  potential 
output  from  these  types  of  vessels,  for  mode  shapes  of  simple  mathematic¬ 
al  form  defined  by  polynomials,  trigonometric  functions,  and  articulated 
elements.  An  ideal  wave  absorber  would  have  a  more  complicated  mode 
shape,  chosen  to  maximize  the  capture  width  at  each  wavelength.  The 
determination  of  this  ideal  mode  shape  has  not  been  made,  but  it  is 
reasonable  to  expect  this  to  be  a  continuous  wavelike  motion  along  the 
body.  The  resulting  power  will  exceed  that  from  any  of  the  simpler  mode 
shapes  illustrated'  in  Section  5.  An  estimate  of  the  ideal  output  can  be 
inferred  from  Figure  15,  which  shows  the  total  capture  width  of  both  the 
trigonometric  and  two-hinge  modeg> reproduced  from  Figures  12  and  13  for 
the  case  of  limited  amplitude  ratio  0.2.  The  envelope  of  these  two 
superposed  curves  is  an  indication  of  the  maximum  output  obtainable  in 
theory  from  an  elongated  flexible  structure.  It  remains  to  be  seen  if* 
a  practical  device  can  be  developed  which  approaches  this  output, 
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Figure  9 :  Properties  of  three  raft  train. 

Curve  (1):  Fraction  of  extracted  energy  at  the  first  hinge. 
Curve  (2)  and  (3):  Vertical  force  at  hinges  1  and  2,  f  = 
F/pgAH. 

Curve  (4)  and  (5):  Relative  rotation  at  hinges  1  and  2, 

| ©2  “  ® 1 1  x  10  and  1 0 ^  ® 2 ^  x  radians. 
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Figure  10:  Efficiency  of  optimum  designs  for  a  double-peaked  spectrum 
Curves  (1)  and  (2):  Three  and  four  raft  trains  for  best 
averaged  efficiency. 

Curves  (3)  and  (4) :  Three  and  four  raft  trains  for  best 
profit . 

Curve  (5):  Spectrum  (not  to  scale). 


